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; ABSTRACT 

We present a multi-wavelength study of galaxy populations in the core of the massive. X-ray luminous cluster XMMU J2235 at 
z=1.39, based on high quality VLT and HST photometry at optical and near-infrared wavelengths. 

We derive luminosity functions in the z, H, and K, bands, approximately corresponding to restframe U, R and z band. These show 
a faint-end slope consistent with being flat, and a characteristic magnitude M* close to passive evolution predictions of M* of local 
massive clusters, with a formation redshift z> 2. 

The color-magnitude and color-mass diagrams show evidence of a tight red sequence (intrinsic scatter < 0.08) of massive galaxies 
already in place, with overall old stellar populations and generally early-type morphology. Beside the red colors, these massive 
(> 6 ■ IO'^Mq) galaxies typically show early-type spectral features, and rest-frame far-UV emission consistent with very low star 
Q^ ■ formation rates (SFR< 0.2Mo/yr). 

' Star forming spectroscopic members, with SFR of up to ~100Mo/yr, are all located at clustercentric distances >250kpc, with the 

, central cluster region already appearing effectively quenched. Most part of the cluster galaxies more massive than 6 • IO'^Mq within 

the studied area do not appear to host significant levels of star formation. 

The high-mass end galaxy populations in the core of this cluster appear to be in a very advanced evolutionary stage, not only in terms 
of formation of the stellar populations, but also of the assembly of the stellar mass. The high-mass end of the galaxy stellar mass 
ky( \ function is essentially already in place. The stellar mass fraction estimated within rsoo (~ 1%, Kroupa IMF) is already similar to that 

. of local massive clusters. 

' On the other hand, surface brightness distribution modeling of the massive red sequence galaxies may suggest that their size is often 

smaller than expected based on the local stellar mass vs size relation. An evolution of the stellar mass vs size relation might imply 
that, in spite of the overall early assembly of these sources, their evolution is not complete, and processes like minor (and likely dry) 
merging might still shape the structural properties of these objects to resemble those of their local counterparts, without substantially 
affecting their stellar mass or host stellar populations. Nonetheless, a definite conclusion on the actual relevance of size evolution for 
the studied early-type sample is precluded by possible systematics and biases. 

Key words, galaxies: clusters: individual: XMMU J2235. 3-2557 - galaxies: evolution - galaxies: luminosity function, mass function 
- galaxies: fundamental parameters - galaxies: high-redshift 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope under Program IDs 10698, 10496, and 10531, and with the 
Very Large Telescope at the ESO Paranal Observatory under Program 



IDs 60.A-9284, 072.A-0706, 073.A-0737, 074.A-0023, 077.A-0177, 
077.A-0110, 079.A-0758, and 081.A-0759 
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1. Introduction 

Galaxy clusters at high redshift provide us with a unique yet 
rare chance to investigate the effect of the highest density envi- 
ronments on the evolution of galaxy populations. Due to the ex- 
treme rarity of massive galaxy clusters, and especially more so at 
early cosmic epochs, even wide-area deep surveys do not probe 
these peculiar environments which thus need to be searched and 
identified in specifically designed surveys, and then followed-up 
with observations across a wide range of wavelengths in order to 
maximize the scientific return of their discovery. 

Thanks to considerable efforts with on-going cluster surveys, 
which are utilizin|; a^variety of methods (e.g., Gladders & Yee 
2000; Rosati et all 1200 2; Eisenhardt et al. 2004; Wils on et al l 
2005; Wittm an et al.ll20 06; Fassbende r.2008fc,Staniszewski et al.l 
2009; Andre on et 311 72009; Demar co et alJ l2010i and refer- 
ences therein), the number of high-redshift clusters is increas- 
ing, however only a handful of spectroscopically confirmed 



galaxy clusters beyond redshift 1 .3 have been discovered to date 



dMullis et al."2005t 'Stanford et al."2 005l l2006t lEisen hardt et all 

l2010l: Tfanaka et aQ 




Iwilsonetal. 2009; Papovich et al 



Galaxy clusters at such high redshifts are not only important 
for studying the emergence of the large scale structure and for 
constraining cosmological models, but also specifically for the 
evolution of massive early-type galaxies which seem to dom- 
inate massive cluster galaxy populations even beyond redshift 
one. The further back in cosmic time we can reach, and thus the 
closer we can get to the major stages of the galaxy formation 
process, the tighter are the constraints we can set on the evolu- 
tion of these galaxies, as well as on the relevance of the cluster 
enviro nment in shaping their physical properties (among 



envrro nment in snapmg tneir pnysicai properties (among many 
others, iToft et all 120041 iBlakeslee et al.l l2006t I AndreonI [20061 



Strazzullo et al. 200( 


\ |De Propris et al. 2007; De Lucia etal. 


2007; Kodamaetal. 


'2007; Holdenetal. 2007; Zirm et aU 



2008; Andreon 2008; Lidman et al. 2008; Gobat et al. 2008; 
Muzzin et al. 2008; Retturaetal. 2010; Mei et al. 2009; 
Colhns et al.. ,2009i; iHilton et ah aOQ'k iRosati et al.i i2009l and 



references therein). 

While the detailed study of galaxies with stellar masses 
M < 10'°Mo is still Umited at z > 1, even with the current 
generation of 10m class telescopes, massive galaxy populations 
can be studied out to high redshift in sufficient detail, in terms of 
their structure, star formation histories, and stellar masses. These 
studies allow us to constrain the evolution of massive galaxies in 
more nuances than the old, and now widely superseded in its 
original form, " monolithic vs hierarchical" question. This pro- 
vides valuable input to theoretical modeling of galaxy evolution, 
the detailed compariso n of cluster and field galaxies being o ne 
of the many examples (iMenci et al.ll2008l:lRomeo et al.ll2008b . 

Taking advantage of the full range of available observations, 
across the widest possible range of cosmic epochs and environ- 
ments, has become particularly important in order to probe dif- 
ferent aspects of galaxy evolution, and to test specific predic- 
tions of theoretical models. As an example, the introduction of 
various forms of so-called feedback mechanisms to modulate 
galaxy evolution has reconciled previously considered "anti- 
hierarchical" observations with hiera rchical predictions (e.g., 
iDe Lucia et al.ll2006tlBower et al.ll2006 . and references therein), 
and elucidated the decoupling between star formation and mass 
assembly histories of massive galaxies. In order to probe these 
two processes independently, galaxy evolution has to be studied 
by directly sampling galaxy populations up to the highest red- 
shifts. 



Fig. 1. A color image of XMMU J2235 obtained by combin- 
ing the ACS F775W and F850LP and NICMOS F160W images 
(North is up. East to the left). The image shows a region of about 
2x2 arcmin^, corresponding to about 1x1 Mpc^ at the clus- 
ter redshift. [Image resolution significantly degraded for astro-ph 
submission] 



In this work we present multi-wavelength observations of 
the X-ray luminous galaxy cluster XMMU J2235-2557 (here- 
after XMMU J2235, RA= 22''35'«20'.82, Dec=-25°57'40.3", 
J2000) at z = 1.39, which is the m ost mas sive virialized struc- 
ture discovered beyond redshift one dJee et a l. 2009; Ros ati et al.l 
2009), spectroscopically confirmed with 30 cluster members 
with redshift in the range 1.37 < z < 1.41. This is a follow- 
up work of the first multi-wavelegth analysis presented in 
Ros ati et al. I (12009), which uses an enhanced data set and is 
aimed at obtaining a more comprehensive picture on the star for- 
mation, stellar mass distribution, and morphological structure of 
the galaxy populations in the cluster central regions. 

Throughout this paper, we adopt a Hq-1Q km s"' Mpc"\ 
Qm=0.3, Qa=0.7 cosmology, and the AB magnitude system. 



2. Data and derived quantities 

2.1. Observations and catalog production 

This work is based on a multi-wavelength dataset collected in 
the field of XMMU J2235 with the VLT and HST in the wave- 
length range 3500A to 2//m. In the following we use VLT 
photometry (VIMOS U, FORS2 R, HAWK-I J and K,) and 
spectroscopy (FORS2), and HST photometry (ACS F775W and 
F850LP - hereafter i and z, and NICMOS F160W - hereafter H). 
A description of the first available HST/ ACS d ata and VLT spec- 
troscopic campaign has been published in IRosati et al.l (120091) : 
in this work we use the full presently avail able ACS dataset, in - 
cluding images from program GO- 10 496 (iDawson et alj|2009h . 
as was presented in Jee et al. (2009"). VLT/HAWK-I data have 
been pubhshed in Lidman et al. (2008), while VLT/FORS2 data 
were presented in Mullis et al. (2005), thus we refer to these pa- 
pers for a full description of these observations. 
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The NICMOS F160W data we used in this work were ob- 
tained in August 2008 (GO 14, Proposal ID 10531, PI: C. 
Mullis), resulting in a mosaic of about 2.5x2.5 arcmin^ with an 
average exposure time of ~ Ihr. 

The U band data (ESO 079.A-0758, PI: M. Nonino) were 
obtained in August 2007 with VIMOS at the VLT. 32 dithered 
observations were collected for a total of 6hrs 20m with seeing 
conditions ranging from 0".5 to V\2. The images wer e reduced 
and stacked in a similar fashion to lNonino et al.l(l2009l) . The zero 
point of the mosaic was derived using standard stars observed in 
photometric nights. The limiting magnitude of the mosaic cov- 
ering the cluster was estimated from the counts distribution in a 
2" apertur43 centered on 5000 random points. After correcting 
for aperture effects (0.19 mag), and for the effect of the correla- 
tion in the noise introduced by the coaddition step (0.2 mags), 
the final value is 28.8 AB mag (Icr). 

Table 1. Summary of the imaging data used in this work. 
Columns 1 and 2 list the instruments and filters with which im- 
ages were acquired. Column 3 gives an estimate of the image 
resolution. Column 4 gives the total area of the image/mosaic 
in each passband, and the actual area used to build the multi- 
wavelength catalog. Further details about these images, and ad- 
ditional information specific to their use for different purposes 
in this work, are given in the text. 



Telescope/Instrument Filter FWHM Area (total/used) 







(arcsec) 


(arcmin-) 


VLT/VIMOS 


U 


0.8 


50/10.3 


VLT/F0RS2 


R_SPECIAL 


0.75 


50/10.3 


HST/ACS 


F775W 


0.1 


11/10.3 


HST/ACS 


F850LP 


0.1 


11/10.3 


VLT/HAWK-I 


J 


0.55 


180/10.3 


HST/NICMOS 


F160W 


0.35 


6.2/6.2 


VLT/HAWK-I 


K. 


0.4 


180/10.3 



Source extraction and photom etry were performed with 
SExtractor ( Berlin & ArnoutsllI99 6'). either in single-image (for 
the determination of MAG_AUTOs used for deriving the lumi- 
nosity functions) or in dual-image mode (for the determination 
of aperture colors used for deriving the broad-band SEDs). We 
adopted MAG_AUTO as an estimate of total magnitude, while 
colors (and thus the broad-band SEDs) were estimated in aper- 
tures of 1" and 1 .5". These aperture sizes where chosen based on 
the typical angular size of galaxies at the redshift we are inter- 
ested in, and on the resolution of the available data. The 1" aper- 
tures allow us to measure colors of cluster members within 1-2 
effective radii (see section|4]i, while for estimating stellar masses 
we used the 1 .5" apertures as a compromise between optimising 
the S/N, minimising the effect of neighbors, and reducing the 
errors on the correction to total masses (see section l23b . 

The angular resolution of our imaging data ranges from a 
FWHM of ~0.1" for the ACS images to ~0.8" for ground-based 
optical images (see Table [TJ. In order to match these different 
resolutions in the multicolor photometric catalog, aperture cor- 
rections were applied in each passband as estimated from the 
growth-curve of point-lik e sources in the fie ld (similar to the 
approach described in e.g. iRettura et al.ir2006l) . This resolution- 

' Here and in the following, aperture sizes refer to the aperture diam- 
eter. 



matching approach allows us to measure accurate photometry 
without degrading the image quality in any of the passbands, 
which is an important advantage when dealing with a crowded 
field such as a cluster corfl 

Mag nitudes w ere corrected for Galactic extinction according 
to Schle gel et al] (U998). 

Multi-wavelength photometric and morphological catalogs 
for the available sample of ~ 30 spectroscopically confirmed 
cluster members will be published in a forthcoming paper 
(Nunez et al., in preparation). 

2.2. Photometric redshifts 

We used the U,R,i,z,J,H,Ks photometric coverage of the XMMU 
J2235 field, to estimate photometric redshifts (photo-zs) by 
comparing the observed photometry with a library of 33 SED 
templates covering a range of star-formation histories, ages 
and dust content. Together with local galaxy templates (e.g., 
Coleman et al. 1980; Mannucci et al. 2001; Kinney et al. 199Q), 
we used a set of semi-empir ical temp l ates based on observations 
plus fitted SED models (iMarastonl Il998t |Bruzual & Chariot 
2003) of ~ 300 galaxies in th e FORS Deep Field (Heidt et aL 
2003) and Hubble Deep Field dWilliams et al.llT9 96). in order to 
better represent objects to higher redshifts. A different prior on 
the redshift distribution is assumed for different types of tem- 
plates (e.g. an old local elliptical template is assumed to be in- 
creasingly unlikely at higher redshifts, while templates corre- 
sponding to young stellar populations or QSOs are assumed to 
have a basically flat likelihood across the whole redshift range 
probed). In addition, a weak broad prior on the absolute optical 
and NIR magnitude lowers the probability to have magnitudes 
brighter than -25 and fainter than -13. Ly-a forest deple tion o f 
galaxy templates is implemented according to Ma daul (fl995l) . 
The "best-fit" photo-z Zphoi is chosen as the redshift maximiz- 
ing the probability among all templates, and an error on Zphot is 
defined as e,pi,ot - [^ijizi - Zphor)^PijV^^, with Zi the considered 
redshift steps, and Pij the contribution of the j-th template to the 
total p r obability function at r e dshift Zj. We refer to Bend er et al.l 
( I2001h ; iGabasch et all (|2004 ; iBrimioulle et"al] (12008) for a de- 
tailed description of both the photo-z estimation method, and 
the construction of the templates. 

Systematic offsets between the measured and template- 
predicted colors as a function of redshift (which may be due for 
instance to errors in the estimated zero-point and aperture cor- 
rections, uncertainties in the filter response curves, or systemat- 
ics in the templates) were estimated using ~ 70 spectroscopic 
redshifts available within the ~ 3' x 3' field, and corrected for 
in order to minimize the systematic shift between observed and 
predicted colors for well-fitted spectroscopic galaxies. 

The photo-z performance against the available spectroscopic 
sample is shown in the top panel of Figure |2] The figure clearly 
shows the impact of the lack of coverage between the U and R 
passbands at lower redshifts, in particular at z<0.6 where the 
4000A break is poorly sampled by the available photometric 



- For aperture sizes used in this work, and for our worse resolution 
of ~ 0.8", we estimate through simulation of synthetic sources (with 
the IRAF task mkobjects) that this resolution-matching approach cor- 
rects for systematics of 0. 1 to 0.5 mag (depending on the aperture size), 
with a residual systematic offset of 0.01-0.04 mag for a de Vaucouleurs 
profile with effective radius 0.24" (~2kpc at z=1.39, about the median 
effective radius of our red-sequence sample), and ~0.05 mag for effec- 
tive radii of 1.2" (~10kpc at z=1.39, about our largest effective radius), 
and a scatter of about 0.01-0.02 mag. 
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coverage, often resulting in unconstrained, inaccurate photo-z's. 
However, by comparison with a spectroscopic sample of ~50 ob- 
jects in the redshift range 1 < z < 2 {~ 40% of these are cluster 
members), we estimate a systematic offset (median Az/(1h-z)) of 
0.006, a scatter of Az/(l-i-z)< 6'%B and 4% catastrophic outliers 
(two objects with |Az/(l-i-z)| > 20%; we note that both sources 
have a spectroscopic redshift deemed uncertain, with a quality 
flag B or worse, see Rosati et al. (2009)). The middle panel of 
Figure |2] shows in more detail the photo-z performance in the 
1 < z < 2 range, suggesting that in the relevant redshift range 
we can derive accurate enough photo-zfl In the bottom panel 
of Figure |2] we compare the spectroscopic (solid black line) and 
photometric (dashed line) redshift distributions in this redshift 
range for bright galaxies (<M*-i-l), where we have enough spec- 
troscopic coverage (see below, also figures in sectionHJ. While 
the peak at the cluster redshift is also visible in the photo-z dis- 
tribution, the gray histogram shows the spread in photo-zs of 
all available spectroscopic members (regardless of their magni- 
tude). 

By comparison with the available spectroscopic sample, we 
estimate that by selecting photo-z candidates within 3cr of the 
cluster redshift (1 < z < 1.8) we include all spectroscopic mem- 
bers, while slightly more than half of the selected sources are 



redshift 
1.0 



interlopers (94% of which in the range 1 < z 



spec 



1.8fl 



2.3. Stellar masses 



Stellar masses were obtained by fitting iBruzual & Charlo3 
(l2003h models to the (resolution-matched) U,R,i,z,J,H, K, SEDs 
as measured in 1 .5" ap ertures (as described in e.g., Go bat et al.l 
l2008HRosati et al.l2 009). Stellar masses are determined for mod- 
els of fixed solar metallicity and with a ISalpeterl (Il955i) IMF, 
however for the purpose of comparison with othe r hterature sam- 
ples (see below) stellar masses are converted to iKroupal (l200l 
IMF masse^. Different determinations of the stellar masses 
were derived, using different star formation histories (exponen- 
tially declining, delayed exponentially declining, with or without 
a recent burst of star formation), including or not dust attenua- 
tion. These different determinations mostly agree within better 
than 50%, (the scatter between different determinations is about 
~30%), and never differ by more than a factor of two. We will 
consider this ~ 30% scatter as an estimate of the typical error 
on the stellar masses due to different SFHs and dust properties, 
however we remind the reader that this estimate does not include 
other systematics, including those due to IMF and metallicity be- 
ing different from what we assumed. 



3 The NMAD estimator jHoaglin et al.ll983l : lllbert et aU2009h gives 
a scatter of 0.045 for the whole spectroscopic sample, and 0.03 in the 
range 1 < z < 2. 

Out of 7 catastrophic failures over the whole < z < 3.5 redshift 
range, four have an uncertain spectroscopic redshift, and two are z~0.3 
sources for which the photo-z is essentially unconstrained, as the error- 
bars in Figure [2] show. 

' While the 3(T photo-z retained samples are used to include all plau- 
sible cluster members, for the purpose of statistical analysis (as for in- 
stance in section |6l(, in the following we also use 2cr photo-z retained 
samples. Based on the comparison with the spectroscopic sample, these 
too are estimated to be virtually 100% complete, and affected by an al- 
most 50% contamination by interlopers. In the stellar mass range rele- 
vant to this work (M> IO'^Mq), the difference between the two samples 
is just ~ 10 sources, with almost 90% of the 3cr sample belonging to 
the 2(T sample as well. 

L0g(MaSSK™»pa) = h0g(MaSSs „l peter) " 0.19. 
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Fig. 2. Top panel: photometric vs spectroscopic redshifts for the 
whole spectroscopic sample within the considered field (one 
point at z=3.26 with photo-z=3.44 is not plotted). The gray line 
traces the bisector. Gray symbols sho w uncertain s pectro scopic 
redshifts (quality flag B or worse, see lRosati et al.l (l2009l) : these 
objects are plotted together with the secure redshifts in the lower 
panels). Middle panel: Az/(l+z) as a function of redshift for the 
spectroscopic sample at 1 < z < 2 (one catastrophic outlier with 
an (uncertain) redshift z= 1 .06 and a photo-z = 1 .64 falls out of the 
plot). The dotted line and the gray shaded area show the median 
and the 16'''-84''' percentile range of Az/(1 H-z). The hatched area 
shows the Az/(1 H- z) scatter as estimated by the NMAD estima- 
tor Bottom panel: the spectroscopic (solid black line) and pho- 
tometric (dashed line) redshift distribution of bright (H< 22.2) 
galaxies in the studied field. The gray histogram shows instead 
the photometric redshift distribution of all spectroscopic cluster 
members. Histograms are slightly offset for clarity. 



In order to obtain total stellar masses, the masses derived 
from the SEDs measured within a radius of 0.75" were renormal- 
ized by the flux ratio between the 1 .5" aperture magnitude and 
the total (MAG_AUTO) magnitude in the z band. This approxi- 
mation neglects the effects of a change in stellar mass-to-light 
ratio in the very external regions of the galaxies. However, given 
the typical galaxy sizes, the correction is generally small: among 
the sample relevant to this work (1 <z< 2, log(M»/M0)>lO.4), 
the median correction is less than 10%, and for the great part 
of the sources is less than 50%. We thus expect this error not 
to significantly affect our results, also in view of the uncer- 
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tainties generally affecting SED-determined stellar masses (e.g., 
iLonghetti & Saraccdl2009h . 

2.4. Structural properties 

We used the GALFIT (iPeng et al.ll2002h software to model sur- 
face brightness profiles of sources brighter than z ~ 24.2. We 
set this magnitude limit based on S/N considerations, according 
to previous results on the accuracy of the retrieved parameters 
as a function of magnitude (S/N) from simulations of surface 
brightness fitting (e.g., among others, Ravindranath et al. 200^ 
iPann ella et al. 2009b, see also below), as well as visual morphol- 
ogy (jPostman et al...2005,) ). This is sufficient to probe the mor- 
phologies of red galaxies in XMMU J2235 down to ~M*h-1. 

We fitted PSF-convolved ISersid (11968 ') profiles to ACS z 
band images, usi ng the PSF mod e ls deri ved from principal com- 
ponent analysis dJee et alJ l2007l l2009l) : using a stack of high 
S/N, unsaturated point-like sources in the field does not change 
our results. The possibility offered by GALFIT to simultane- 
ously fit multiple sources is particularly useful in a cluster en- 
vironment: when fitting a galaxy surface brightness, all nearby 
source^ brighter than magg^j^j^^y-i-S were also modeled at the 
same time, reducing biases produced by contamination by neigh- 
bors. 

A detailed description of the morphological properties of 
bright cluster galaxies, and a comprehensive morphological 
analysis including the high-resolution ground-based NIR data, 
will be presented in Nunez et al. (in preparation). In the follow- 
ing we only use the Sersic index tisersic to broadly classify galax- 
ies in early and late morphological types, and the estimated (cir- 
cularized) effective radius as an estimate of the galaxy size. 

3. Luminosity functions 

The luminosity function (LF) of galaxies in the central region 
of XMMU J2235 was determined in three passbands (z, H, K^), 
roughly corresponding to restframe U, R, and z bands. The re- 
gions actually used for the LF determination were about 7.5, 6, 
and 7 arcmin^, in the z, H, and K, bands, respectively. At the 
cluster redshift, the area probed extends out to a clustercentric 
distance of about 700 kpc. The completeness limits of the z, H, 
and Kj band images was estimated based on the turnover of the 
number counts of sources with S /N > 10, and turn out to be ~ 
25.3, 25, and 23, respectively. 

For the purpose of statistical subtraction of background con- 
tamination, we also made use of publicly available photometry 
acquired with the same, or very similar, instrument/filter as the 
XMMU J2235 images. For the band we used VLT/ISAAC 
Kj band photometry in a ~ 60 square arcmin region within 
the GOODS CDF-S field (v2.0 release, Retzlaff et al. (2010)). 
Only regions with a lOcr depth comparable or better than the 
XMMU J2235 Kj image were used. For the H band we used the 
NICMOS im aging on a ~ 5arcmin ^ region in the Hubble Ultra 
Deep Field (Thom pson et al .'2005), while for the z band we used 
ACS ph otometry in the G OODS (North and South) fields (v2.0 
release, iGiavalisco et al.l (|2004), Giavalisco et al. in prep.), for 
an overall area of more than 200arcmin^. 

For all fields, point-like sources were identified and removed 
based on the combination of morphological/concentration 
parameters estimated by SExtractor (MAG_AUTO, 

' Sources within a region of dimensions about 4 times the dimension 
of the source as estimated through SExtractor parameters (ISO_AREA, 
position angle). 



FLUX_RADIUS, CLASS_STAR, FWHM.IMAGE) on the 
z band image. 

Luminosity functions (LFs) were determined by means of 
statistical background subtraction, but at the same time taking 
full advantage of the available spectroscopic and photometric 
redshift information. The best purely statistical approach would 
avoid binning the data, and would derive the cluster galaxy LF 
by simultaneously modeling the number counts in the cluster 
and in a control field, as described in detail in lAndreon et al.l 
( 12005 *). However, this approach would not allow us to take into 
account the information about cluster membership that we derive 
by other means (spectroscopic and photometric redshifts). While 
information on cluster membership is obviously always relevant, 
its specific importance in the determination of the LF (and its er- 
rors) depends on different factors including the quality of such 
membership information and how much it contributes to the con- 
straints on the LF. When only one cluster field is available, the 
statistics at the LF bright end is very poor; furthermore, partic- 
ularly in the case of very distant clusters, the imaging depth can 
reach at most about 3-4 mag fainter than M*. In these conditions, 
it is very important to make full use of all the available informa- 
tion, in order to best constrain the LF even when the intrinsically 
low counts would make its purely statistical determination rela- 
tively loose. 

We thus adopted the following approach (wh ich extends with 
the us e of photo-zs the approach adopted in IStrazzullo et al.l 
(l2006l) ): we first determined the LF in each passband by sub- 
tracting the (area normalised) counts in the control field from the 
counts in the cluster field. The eiTor on the excess counts was de- 
termined by summing in quadrature the Poissonian errorfl from 
both cluster and field counts (even though it is dominated by 
the cluster contribution when the cluster field is much smaller 
than the control field). The eiTor contribution coming from cos- 
mic variance due to galaxy clustering, as estimated according to 
lHuangetal.l(ll997h . is much smaller than the Poisson error (a 
factor of a few percent at most) and w as neglected (see also e.g., 
lAndreon et al.ll2005l: IStrazzullo et al.l l2006). 

We then took into account the spectroscopic information, im- 
posing that in each magnitude bin the excess counts are at least 
equal to the spectroscopic members in that bin. The spectro- 
scopic (and photo-z) information is also considered in the deter- 
mination of the final error on the excess counts in each bin. The 
membership for the brightest galaxy populations in the core of 
XMMU J2235 can be considered quite well established (see also 
figures in section |4]below). In the critical (for the LF determina- 
tion) range from the BCG magnitude to M*, the spectroscopic 
completeness is about 70%. Furthermore, basically all plausi- 
ble members brighter than M* were targeted for spectroscopic 
follow-up, and thus the remaining sources without spectroscopic 
redshift have a photo-z far below the cluster redshift, and well 
beyond 3cr from the cluster redshift according to the photo-z 
estimated rms. Therefore, in estimating the error on the excess 
counts, when all objects in a magnitude bin have a spectroscopic 
redshift, or a photo-z beyond 3cr from the cluster redshift, we 
assumed that: i) the LF is basically "membership determined" 
without the error contribution from the background counts, ii) 
the lower limit on the excess counts cannot go below the lower 
limit based on the spectroscopic excess counts, and iii) the up- 
per limit on the excess counts is limited by the upper limit on the 
counts of spectroscopic members plus photo-z candidate mem- 
bers within 3cr from the cluster redshift. We remind the reader 



Following [Geiirel^ ( fToM) , upper limits are estimated as A' + 1 + 
(N + 0.75)"^ and lower limits as N(l - ^ - 
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Fig. 3. Bottom panels: The luminosity function of galaxies in the central region (within ~ 700kpc of the cluster center) of XMMU 
J2235, in the z, H and K, bands. Solid symbols show binned counts (with Icr errors, see text for details), and black lines show 
the best-fit Schechter function as determined for magnitudes brighter than ~M*+2. Confidence levels (1,2, 3cr) on the two relevant 
parameters (characteristic magnitude M* and faint end slope a) are plotted in the small insets. In the central panel (H band LF) 
and corresponding inset, the gray lines show the best-fit Schechter function and confidence levels for the fit down to ~M*+4. 
The formal best-fit values and Icr errors for M* and a, as fitted down to ~M*+2, in the three passbands are: M*,=22.8 + 0.6, 
a,=-0.95 + 0.6, M*H=21.4+i 5 , aH=-l!i.4 ( MV=21.2+JJ;5 , aH^-l.2^_°A5 down to ~M*+4), MV,^20.9^Q;^, ay.= -l>!;i Upper 
panels: The evolution of M* with redshift as observed in the z, H and band, according to Kodama & Arimotol (1199 7) models, 
for three formation redshifts Zf-2,3,5. Solid symbols show M* (and errors) as determined from the LFs in the lower panels. 



that the 3cr photo-z retained sample is estimated to be 100% 
complete, and affected by a >50% contamination, thus we ex- 
pect the upper limits based on such sample to be robust. 

In Figure [3] we show the binned LFs determined in this way 
in the three z, H, and passbands. The available photometry 
reaches out to ~M*-i-2 or M*+4 depending on the passband. For 
ease of comparison, all LFs are plotted in Figure [3] over a mag- 
nitude range of ab out 5 magnitudes , from ~M*-2 to ~M*-i-3.5. 
The formal best-fit lSchechteii d 19761) function determined by 
minimization on the binned counts plotted, with errors as de- 
scribed above, is also shown. For the H band we plot both the 
LF determined for galaxies brighter than M*-i-2, as for the z and 
Kg bands, and down to ~M*-i-4, which is only probed by the H 
band. In all cases, the best-fit faint-end slope is close to flat and 
the characteristic magnitude M* is close to the predictions of 
passive evolution of the local M*, assuming that the bulk of the 
stellar populations are formed at a redshift ~ 3. This is shown in 
the upper pa nels of Figure [3] sho wing the redshift evolution of 
M* based on lKodama & Arimotol (119971) models, in each of the 
three passbands, together with the M* determination from the 
LFs in the lower panels. The errors on M* are determined based 
on the two-parameter confidence levels plotted in the small in- 
sets in the lower panels. 

While the band still samples the restframe near-infrared 
light, which can be considered as a probe of the stellar mass, the 



H and particularly z bands sample wavelengths more affected 
by recent or on-going star formation activity. Nonetheless, 
the LFs in all three passbands suggest a similar M* evolu- 
tion. The flat faint-end slope, and the measured M* close to 
passive evolution predictions of the local cluster galaxy LF, 
point toward the LF bright end being similar to that of lo- 
cal clusters, beside evolution of the stellar populations. This 
extends to z ~ L4 previous findings on the early assembly 
of the massive galaxy populations in the cluster core regions 
(DeProprisetal. 1999; Nakataetal. 2001; Kodama & Bow3 
200 31: iToft et all 12 003; ElHs & Jones 2004; Toftetal. 200l 
AndreonI 120061; ' Strazzullo et al., ,2006(: i De Propris et ah j2007l : 
Muzzin et al, 2008; Manco ne et al.ll20ld)~ 



4. Stellar population properties of cluster galaxies 

In the following, we investigate the properties of cluster galaxies 
in XMMU J2235, and of their host stellar populations, by means 
of direct, simpler approaches (as the color-magnitude diagrams) 
and of more detailed modeling (SED fitting with stellar popula- 
tion synthesis models), using the photometric and spectroscopic 
measurements described above. 
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Fig. 4. The color-magnitude diagram in the central region (r < 
600kpc) of XMMU J2235. All sources in the field are shown, 
except for point-like sources. Filled circles are spectroscopic 
members and large empty circles show photo-z retained sources 
(0.96 < z < 1.82). Gray symbols (crosses and pluses) show 
spectroscopic and photo-z interlopers, respectively. The colored 
lines show th e expected locat i on of the red sequence at z ~ 1 .4 
based on Ko dama & Arimot ol dl997h models for different forma- 
tion redshifts (z/=2, 3, 5, as indicated). The slanted continuous 
line shows the IQcr completeness (for the overall catalog, con- 
sidering AUTO apertures). The plotted z-H colors are derived 
from resolution-matched 1" aperture magnitudes, and error bars 
include the estimated error contribution from the resolution- 
matching procedure. Error bars are plotted only for retained 
sources. 



4.1. Color-magnitude diagrams 

Color-magnitude diagrams (CMDs) shown in the following re- 
fer to a region of about 3 x 3arcmin^ in the cluster center (the 
maximum distance from the BCG is about 1.2Mpc). The H band 
image only covers a portion of this area, corresponding to a 
~ 2.4 X 2.4arcmin- square centered on the BCG, which trans- 
lates into ~ 1200 X 1200kpc2 at z=1.39. 

Colors used in the following are based on resolution- 
matched magnitudes measured in a 1" aperture, which minimize 
the photometric errors (and possible contamination by neigh- 
bors), especially for fainter sources. We note that an indepen- 
dent, alternative approach based on model total magnitudes and 
aperture colors obtained via surface brightness fitting of each 
source, gives results in very good agreement with those obtained 
here (Nunez et al., in preparation). We thus assume that the 
systematic errors affecting MAG_AUTO as an estimate of total 
magnitude, and aperture colors obtained as described above, do 
not significantly affect our results. 

In Figure m we plot the z-H vs H CMD. All objects in the 
field are plotted, except for point-like sources which were re- 
moved (see section O. Spectroscopic cluster members are plot- 
ted as filled symbols, while spectros copic interlopers are plot- 
ted as cross e s. For consistency with iLidman et al.l (l2008h and 
iRosati et"al] (l2009l) . we define as cluster members all objects 
with 1.37 < z < 1.41. 

In addition to excluding spectroscopic interlopers, we can 
further clean the CMD in Figure |4] of obvious fore- and back- 
ground objects excluding sources with a photo-z beyond 3cr of 
the cluster redshift (photo-z< 0.96 or > 1.82). These most likely 



interlopers are plotted in Figure |4] as plus symbols, while the re- 
tained sources with photo-z within 3cr of the cluster redshift are 
plotted as empty circles. The whole H<24.5 galaxy sample plot- 
ted in Figure |4] contains 391 sources, out of which 205 have a 
photo-z within 3cr from the cluster redshift, 21 are spectroscopic 
cluster members, and 24 are spectroscopic interlopers. As ex- 
plained above, the sample selected by retaining photo-z candi- 
date members within 3cr of the cluster redshift is virtually 100% 
complete, but affected by a > 50% contamination, thus we point 
out that the only purpose of the "photo-z cleaned samples" plot- 
ted in the following Figures |5] and |6] is to show more clearly the 
bright cluster population by removing the contamination of ob- 
vious foreground sources. 

In Figure|5] we plot the z-H vs H and z-J vs J CMDs, show- 
ing only spectroscopic cluster members and photo-z retained 
sources. The observed z band well matches the restframe U band 
of the cluster galaxies, and the J and H bands sample the rest- 
frame SED of cluster members at ~5000A and ~6700A, respec- 
tively. The colors shown in Figure |5] thus correspond approx- 
imately to restframe colors U-B/U-V and U-R. While the ob- 
served z-J color best samples the amplitude of the 4000A break, 
the z-H vs H CMD benefits from the excellent photometric accu- 
racy and depth provided by HST imaging. Nonetheless, the pres- 
ence of a tight red sequence is evident in both CMDs, with the 
zero-point and slope of the bright red sequence in good agree- 
ment with iKodama & Arimotor(ll997h model predictions for a 
formation redshift Zf ~ 3. The intrinsic scatter of the bright z-H 
red sequence (within the shaded area of Figure |5]l is estimated 
to be 0.08+0.01, independent of considering only morphologi- 
cal early-types or all red sequence galaxies. This estimate of the 
intrinsic scatter only takes into account photometric errors as es- 
timated by SExtractor; considering the error introduced by the 
adopted PSF matching approach could lower the intrinsic scatter 
to values as low as 0.05±0.015. 

We can roughly estimate the backward evolution of the 
XMMU J2235 red sequence bright-end (within the shaded area 
in Figure|5]) with a simplistic approach. The restframe U-B scat- 
ter calculated from synthetic colors for the red-sequence galaxies 
is ~ 0.06, in goo d agreement with previous es timates a t simi- 
lar redshifts (Blak eslee et"al]|2006l: lOobat et al. l 2008; M ei et all 
2009). Assuming that the star formation histories of the bright 
red sequence galaxies may be described by a simple exponen- 
tially declining star formation rate with time-scale t, we evolved 
the z ~ 1.39 red-sequence galaxies to earlier epochs (see also a 
more detailed analysis in Gobat et al. (2008)). In general agree- 
ment with the IGobat et alj (12008 ) results based on the red se- 
quence of RDCS1252 at z ~ 1.2, we find that beyond redshift 2 
(e.g. ~lGyr earlier than the observed epoch for XMMU J2235) 
a significant part of the red sequence galaxies would no longer 
be classified as "red", and the scatter of the remaining "red se- 
quence" would be substantially larger than what measured in 
XMMU J2235. We note that this simplistic de-evolution of the 
red sequence does not take into account the whole population of 
massive galaxies in XMMU J2235 and the effects of dust extinc- 
tion, thus does not include the contribution to the red sequence 
by reddened starbursts. By redshift ~2.3, ~ 40% of the initial 
red-sequence sources would not be "red"; those still appearing 
as a "red sequence" would have an average age/r ~5, thus they 
would have joined the red sequence relatively recently. This is 
in overall agreement with red-sequence ob servations at z > 2 in 
proto-clusters as well as in t he fiel d (e.g., iKodama et al1l2007l; 
iKriek et aUlIOOStlZirm et al.ll2008h . 
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We note that these CMDs are based on resolution-matched, 
fixed-aperture photometry in a 1" aperture, corresponding to 
a physical radius of ~ 4.2kpc, and to about 1.8 times the 
median effective radius measured for the bright red-sequence 
early-types. Since the resolution-matching is performed inde- 
pendently on each passband image, the agreement of the forma- 
tion epochs as estimated from the z-H and z-J CMDs suggests 
that the determined CMDs should not be significantly affected 
by systematics introduced by the resolution-matching proce- 
dure. 

The presence of a ti ght red sequence in X MMU J2235 has al- 
ready been reported in iLidman et al.l ("2008), based on the same 
high quality ground-base d NIR photometry used in this work. 
The iLidman etaP (l2008l) J -Ks vs Kj CMD was built within a 
distance of 900 kpc from the cluster center, and the parameters 
of the red sequence were determined based on 9 bright red- 
sequence galaxies within a smaller region (r < 90kpc from 
the cluster center), four of which are spectroscopic members. 
Zero-point and slope of the J-K^ vs red sequence are in 
good agreement with simple stellar population (SSP) models, 
reproducing the red sequence in the Coma cluster, formed at 
z ~ 4. Therefore, while the observed J-K, color at z=1.39 is 
not the best choice to probe the stellar population ages, and 
while even high quality ground-based imaging hardly competes 
with HST imaging, especially in crowded fields, our z-H vs H 
color-magnitude rel ation substantially confirms the results of 
ILidman et al.l (l2008l) on the presence of a tight red sequence of 
bright, passively evolving galaxies which formed their stars at 
high redshift. In addition, we note that the formation epoch of 
massive red sequence galaxies estimated in Figure |5] by com- 
parison with Kodama & Arimoto models, are in excellent agree- 
ment with the star-formation history derive d by modehng the 
spectro- photomeric data of th ese galaxies jRosati et al.l l2009l) 
with Bru zual & Chariot! (l2003h models. 

In Figure |5] spectroscopic cluster members are color-coded 
according to the detection or non-detection of [Oil] emission 
in their spectra; however, note that due to the quality of our 
spectraJOII] equivalent widths of less than 5A cannot be de- 
tected dRosa ti et al. 2009). As noted in Lidman et al. (2008) and 
iRosati et aP ^2009), emission line members tend to avoid the 
central area of the cluster, the closest emission-line member be- 
ing at a distance of ~ 240kpc from the BCG. We note that, even 
though the spectroscopic sample is generally not complete, ac- 
cording to the photo-z selection it is likely to be complete (in 
the central ~ 1200 x 1200kpc^ region) for sources brighter than 
M*, and > 80% complete for masses larger than 1O"M0). 

We also make use of the morphological information in the 
ACS z band to broadly classify sources in early and late mor- 
phological types. We are mainly interested in identifying late 
type galaxies on the red sequence, since these might be dusty 
star-forming galaxies landing on the red sequence because of 
dust reddening instead of the old age of their stars. As a re- 
sult, these objects might bias our analysis and interpretation of 
the red sequence. The separation of early and late morpholog- 
ical types based on the Se rsic index riser sic alon e is known not 
to be very accurate (e.g., iBlakeslee et al.ll2b03h . and a multi- 
parametric classification should be adopted for a robust char- 
acterization. However, in this work we only aim at classifying 
galaxies in very broad morphological classes, which can be ac- 
complished with g ood statistical accuracy by using the Sersic 
index alone (e.g. Blant on et aLll2003l: iRavindranath et al.ll2004t 
[Pannella et al. 2006). In the following, we will use risersic > 2 
as the threshold to identify morphological early-types. With this 
selection, we are likely to include basically all ellipticals, and 
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Fig. 5. Color-magnitude diagrams in the central region of 
XMMU J2235. Only spectroscopic members and photo-z re- 
tained sources (0.96 < z < 1.82) are shown. Filled blue and 
red symbols show spectroscopically confirmed members with 
and without detectable [Oil] emission, respectively. The dotted, 
dashed and long-dashe d lines show the expe c ted lo cation of the 
red sequence based on iKodama & Arimot ol ([T997h models for 
different formation redshifts. The dark purple line in the lower 
panel shows a fit to the red sequence galaxies within the gray 
shaded region. The shape of all symbols are coded according to 
morphological classification in the z band: circles are galaxies 
with rtsersic > 2, spiral symbols are galaxies with tisersw < 2, 
squares are (mostly faint) objects for which no reliable classifi- 
cation is available. 



the vast majority of SO galaxies, while having some contamina- 
tion mainly from late-types with a significant bulge component. 
Such a threshold is a conservative choice from the point of view 
of the characterization of the red sequence sample, which will 
include some late morphological types. On the other hand, when 
considering the morphological mix of galaxy populations on the 
red sequence, we remind the reader that our early-type sample 
might be affected by such contamination. 

In Figure|5] we thus plot with different symbols galaxies for 
which we were able to perform a morphological analysis, ac- 
cording to their Sersic index (circles for nsersw > 2, spiral sym- 
bols for nsersic < 2). Most of the bright red-sequence galaxies 
in the cluster core are morphologically early types according to 
this criterion (~ 90% down to M*). Visual inspection of the im- 
ages of red sequence galaxies (gray shaded area in Figure |5]) 
confirmed the results based on the automated classification, in 
particular a dominance of early-type galaxies and a minor con- 
tribution from late types (often hosting a prominent bulge). 

All but one of the spectroscopic members classified as mor- 
phological early types do not have detectable [Oil] emission. 
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The one exception is a galaxy with risersic ~2.4 and intermedi- 
ate colors (z-J~l.l,z-K5 ~1.8), located more than 900kpc away 
from the cluster center. [Oil] emission is detected in most of the 
spectroscopic members classified as morphological late types. A 
fraction (~ 30%) of the (incomplete) sample of spectroscopic 
members classified as morphological late types lie on or close to 
the red sequence. Also, about 40% of the spectroscopic members 
morphologically classified as late-types within the studied field 
have no detectable [Oil] emission, thus likely an [Oil] EW<5A. 
All but one of these late-type galaxies without detectable [Oil] 
are relatively blue (as compared to the red sequence); none of 
them is located in the very central area of the cluster, the closest 
to the BCG is at a distance of about 150 kpc, while the others 
are more than 600 kpc away. 

Unfortunately, our sample of late-type spectroscopic mem- 
bers is inevitably incomplete (many late-type candidate mem- 
bers are beyond the spectroscopic limit). Much larger sam- 
ples will be needed to disentangle morphological and spectro- 
photometric evolution of cluster galaxies at this redshift. We also 
note that [Oil] emission might also be due to AGN activity rather 
than star formation ( Yan et al. 2006) , especially for objects with 
a typically red SED. Nonetheless, we also note that the only 
[Oil] emitter actually lying on the red sequence has a late-type 
morphology; its colors would be ~ 0.2 - 0.3mag bluer in a larger 
aperture including the extended disk. This source is located at a 
distance of ~ 860kpc from the cluster center. None of the red- 
sequence galaxies shown in Figure |5] is detected as an X-ray 
point source. We note that the 190 ks Chandra observations of 
this field allow the detection of AGN with X-ray luminosities 
Lx > 10^^ erg/s in the 2-10 keV band. 

4.2. Color-mass diagrams and star formation rates 

We show in Figure|6]the color-stellar mass diagram in the central 
cluster area, as in Figure |5] While overall similar to the color- 
magnitude diagram in terms of bimodality of the galaxy popu- 
lations, the color-mass diagram yields more fundamental insight 
because of the disentanglement between color, magnitude, and 
mass-to-light ratio. 

Colors were derived from resolution-matched photometry 
measured in 1" apertures, as described above. The derivation 
of stellar masses is described in section 12.31 and an approxi- 
mate estimate of the mass completeness limit of this sample, 
derived from the z band lOcr completeness assuming a SSP of 
solar metallicity formed at 5 < z < 10, is shown by the thick 
gray line in Figure |6] 

Figure|6]confirms that high mass galaxies (> IO^'Mq) are al- 
ready in place in the core region of XMMU J2235, and that this 
high-mass population is dominated by red-sequence galaxies, 
mainly of early-type morphology, generally lacking evidence 
of residual star formation activity. The stellar mass of the BCG 
is cl ose to th e stellar mass of local BCGs in massive clusters 
(see lColUnset a l. 2009; Stott et al. 2010, for a broader picture) . 
According to semi-analytical predictions dPe Lucia et al ] l2006h . 
a BCG at z~1.4 will have typically assembled only ~ 20% of 
the final stellar mass it will reach at z=0. In such a model, the 
BCG of XMMU J2235 would thus eventually become an object 
of 2.3-lO'^M0; while not impossible, this would certainly put it 
at the very high-mass end of local BCGs. We note however that 
this same model predicts a slower evolution for BCGs of most 
massive. X-ray luminous clusters as XMMU J2235. 

Figure |6] clearly shows that, even at z ~ 1 .4, there are es- 
sentially no blue galaxies in the most massive populations liv- 
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Fig. 6. Color-mass diagrams in the central region of XMMU 
J2235. Only spectroscopic members and photo-z retained 
sources (0.96 < z < 1.82) are shown. Symbols and color coding 
as in Figure |5] The thick gray line across all panels shows the 
estimated mass completeness limit (see text). 



ing in the core of such a massive cluster As expected, going to 
lower masses the importance of late-type galaxies (from both 
the morphological and the photometric point of view) increases. 
It may be interesting to note that blue, disk-dominated spectro- 
scopic members with [Oil] emission generally lie at masses be- 
low ~ 5 X 10'*'Mq. We note again here that the spectroscopic 
sample is biased against red galaxies at these masses, because of 
their intrinsic faintness due to their higher mass-to-light ratio, 
and to the lack of emission lines. 

While keeping in mind that the rest-frame FUV light may 
be heavily affected by dust attenuation, we further note here 
that the rest-frame 1500A flux of massive (> 6 x IO'^Mq) 
red sequence early-types is indeed consistent with these galax- 
ies hosting no significant star formation. After removing galax- 
ies whose U band photometry is likely contaminated by other 
sources, we stacked U band images for eight early-types more 
massive than 6x lO'^'M© (six out of these eight are spectroscopic 
members). We found no significant detection and obtained a 5cr 
U-band upper limit of 27.5 AB mag, corresponding to a (not 
dust-correct ed) star formation r ate (SFR) of ~O.15M0/yr (for 
comparison, ' Rettura et al.l (1201 Ol) estimated by U-band stacking 
an upper limit of ~O.3M0/yr for the SFR in massive early-types 
in the z ~ 1.2 cluster RDCS J1252-2927). 

For comparison, the SFR of the spectroscopic members lying 
below the Zf -2 model in Figure |5] was also estimated based on 
the measured U band (rest-frame ~ 1500 A) photometry. Since 
we know that these are star-forming sources, we estimated a dust 
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attenuation correction for the ~ 1500A res t-frame luminosity 
from the restframe far-UV - near-UV color (" Paddi et al.ll2007t 
[Salimetal. 2007; Pannella et al. 2009a), which in our case in 
sampled by the observed U and R bands. SFRs derived with 
these dust attenuations following e.g. Daddi et^ (l2007h . range 
between 6 and SO-QOMo/yr depending on the assumed relation, 
with a median of ~2OM0/yr. Specific SFRs (SSFR, SFR/stellar 
mass) are between 0.2 and 2-6Gyr"' (also depending on the spe- 
cific relation used), with an average of ~lGyr We remind 
the reader of the uncertainties which affect these values due 
to the scatter in the assumed relations, filter mis-match, and 
general difficulties in estimating coiTec t dust attenu ation values 
for a given source/sample (for instance. ISalim et aP (2007) note 
that their attenuations as a function of U V color are systemat - 
ically lower than other estimates by e.g. iMeurer et al.l ( Il999h : 
ISeibert et alJ (2005)). For comparison we also estimated dust ex- 
tinction for these sources with a different approach, based on the 
correlation between dust attenuation and stellar mass for star- 
forming galaxies, assuming that these objects follow a correla- 
tion similar to that of field galaxies at comparable redshift. We 
estimated the dust attenuation for the 1500A rest-frame lumi- 
nosity based on the galaxy stellar mas s according to t he rela - 
tion determined for z ~1.7 galaxies in iPannella et al] ( l2009ah . 
obtaining an independent estimate of the SFR. We note that 
since the relation between stellar mass and dust attenuation is 
redshift dependent (given the secular decline of (S)SFR), we 
corrected for the cosmic tim e diff'erence be t ween z =1.7 and 
z = 1.39 according to eq.2 in iPannella eTaP (l2009ah . The SFRs 
for the blue cluster members derived in this way range between 
~ 6M0/yr and ~ lOOMo/yr, with a median of ~ 25M0/yr. 
SSFRs range between ~0.3Gyr ' and ~ 3Gyr with an aver- 
age SSFR~1.3Gyr '. Stellar masses of these sources are in the 
range 1-7-10'"Mq. Due to the observational bias, most of the 
passive spectroscopic members lie at masses above this range. 
At masses below ~ 5 ■ lO'^M©, with the available spectroscopic 
data, redshifts can be measured only for sources with a relevant 
amount of star formation (and without extreme dust extinction). 
As Figure|6]shows, the star forming spectroscopic members dis- 
cussed above are among the bluest sources in the cluster galaxy 
populations in their mass range. 



Finally, Figure Q summarizes the information concerning 
stellar masses, (spectro-)photometric and morphological prop- 
erties of 2cr photo-z retained candidate cluster members to- 
gether with their projected distance from the cluster centre. 
Recent work on high redshift clusters at z ~ 1.45 and 1.62, 
corresponding to cosmic times about 0.2Gyr and 0.6Gyr ear- 
lier than XMMU J2235, suggested that a significant star for- 
mation activity mi ght be occurring in the very centra l regions 
(r < 250 - 500kpc. lHilton et alJl2010HTran et al.ll20"Toh of clus- 
ters in a less advanced evolutionary stage. However, as Figure 
|7] shows, star formation appears to be effectively quenched in 
the core of a more evolved, massive X-ray luminous cluster as 
XMMU J2235. The relevance of very reddened star formation 
contaminating the cluster red sequence could be addressed with 
future observations (NIR spectroscopy, Herschel IR photome- 
try), however, even though the rest-frame NUV (U band) pho- 
tometry might be affected by dust extinction, the early-type mor- 
phological appearance and spectral features of the central mas- 
sive red galaxies suggest that these are more hkely passive early- 
types than highly attenuated starbursts. 
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Fig. 7. The projected distribution around the cluster BCG of 2-cr 
retained photo-z candidate members (spectroscopic interlopers 
are not plotted). Symbol shape is coded as in Figures|5]and|6] ex- 
ternal color coding reflect the z-H color while symbol size scales 
with stellar mass (as indicated). Color filling, for spectroscopic 
members, is blue for galaxies with detected [Oil] emission, red 
otherwise. The small inset at the bottom of the plot shows a 
larger view of the central ~ 0.2 x 0.2arcmin^ (~ 98 x 98kpc^). 
Gray lines show radii of 100 to 700 kpc (with a step of 100 kpc). 



5. Stellar masses and sizes of massive cluster 
early-types 

The above analysis in the central regions of a particularly mas- 
sive, distant X-ray luminous cluster shows that the high-mass 
population of cluster galaxies primarily of early-type morphol- 
ogy is substantially in place at epochs earlier than redshift one, in 
relation to both the formation of the bulk of their stars (at z > 2) 
and of the assembly of the bulk of their stellar mass. This points 
toward little evolution occurring in massive galaxy populations 
in the most dense environments over the last 9 billion years. On 
the other hand, it is interesting to remember that the structural 
properties of such massive objects might be different from those 
of similar galaxies in the local Universe. As pointed out in sev- 
eral studi es (e.g., Daddi et al. 2005; Truiiflo et al. 2006b a, 2007; 
Longhetti et al. 2007; Cimatti et al. 2008; Buitrago et al. 20o5 
van der Wei et al.. ,2C)08: .Damianov et al. .20()9; Rettura et al] 
20101: IWilliams et alJlIolOt iTavlor et alJl2010l and references 
therein), massive early-type galaxies at high redshift appear 
to be (on average) more compact than their local counter- 
parts, their size being smaller than that of local early-types of 
compara ble mass (b u t see also resu l ts and discussion in e.g. 
Mancini et alj (1201 Oh ; lOnodera et"al] (l2010h ; IValentinuzzi et al.l 
(.201Q) ;iSaracco et alJ J2OIOI) . and references therein) Indeed, as 
shown in Figure [8] also the massive early-type galaxies in the 
core of XMMU J2235 appear to be generally smaller than those 
in the local Universe. 
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Figure[8]shows the stellar mass against the (circularized) ef- 
fective radiufl as measured from the surface brightness fit in 
the z band (section 12.41 ). for bright early-typ43 red sequence 
galaxies (within the gray shaded area in Figure |5]). Most of this 
sample is made of spectroscopic cluster members. For compari- 
son we show the local stellar mass vs size rel a tion as measured 
for SDSS eai-ly-type galaxies by IShen et alJ (12003 5. which is 
commo nly used as the lo cal reference. We note that the size 
used by IShen et aL I (f200l was measured in the z band, while 
in our case galaxy sizes are measured in the observed z band, 
corresponding to the restframe U band, which might imply is- 
sues due to morphological K-corrections. However, sizes mea- 
sured in bluer bands are ex pected to be larger than sizes mea- 
sured in redder ban ds (e.g ., iBarden et al.ll2005t iMcIntosh et al.l 
120051: iTruiillo et al.ll200 7). Therefore our galaxy sizes, if mea- 
sured in the restframe z band would be, if anything, smaller than 
those plotted in Figure [8] as it is indeed suggested for this very 
sample by the comparison of surface brightness fitting results in 
the observed z and Kj passbands (Nunez et al., in preparation). 
Therefore, taken at face value, the comparison of our early-type 
sample with the local Shen et al. (2003) reference implies that 
the size of massive early-types in XMMU J2235 is on average 
about ~50% of what expected based on the local relation (me- 
dian r/r-=o=0.46±0.08). We note that the most massive object in 
Figure[8]is the cluster BCG; while the size of this galaxy is quite 
large and it lies on the local relation, we note that BCGs tend 
to be larger than similarly massive (non-BCG) galaxies (e.g., 
I von der Linden et al.ll2007l) . 

The smaller sizes of XMMU J2235 massive early-types 
compared to local early-types suggest a scenario where the evo- 
lution of these massive red-sequence galaxies is not actually 
fully completed at early epochs, as would be suggested by the 
color-magnitude and color-mass relations, as well as the LF 
evolution. This leaves room for processes such as minor (and 
likely fairly dry) merging events to be relevant at later epochs 
to shape the final structural properties of these galaxies, without 
substanti ally altering th e ir stel lar masses and stellar populations 
(see e.g. 'Hopkins et al.' ('2010') and references therein; see also 
discussion in Nipoti et al. (2009)). 

On the other hand, we should note that the comparison of 
the stellar mass vs size relation with a local sample, as shown 
in Figure [8] may be affected by systematics in the determina- 
tion of stellar masses at different redshifts, observational bi- 
ases in the determin ation of galaxy sizes (Pannella et al. 2009b; 
iMancini et al. I I20T0I) . and a possible mismatch between the ob- 
served high redshift sample and the local reference, which are 
discussed in more detail below. 

In Figure |8] we also plot the local stellar mass vs. size 
relation for a sample of ne arby cluster galaxies as measured 
by [Valentinuzzi et al.l (I2010h . While consistent at Icr with the 
IShen et al.i (|2003D relation, the IValentinuzzi et al.1 (1201 Ol) de- 



' Six (one) sources out of the 15 plotted in Figure [8] have a for- 
mal uncertainty on their size lower than 15% (10%). Keeping into 
account results from simulations or multiple independent fittin g (e.g. 
iPignatelli et al.l2006l:lBu"itrago et alj2008l ; lPannella et al.l2009bi) on the 
typical uncertainty on the galaxy size (~ 10-15%) achievable by surface 
brightness fitting for sources of S/N similar to those used here, we in- 
crease the error bar for these six sources to reflect an uncertainty of 
15%. 

Galaxies in this sample were classified as "early-types" based on 
their Sersic index larger than 2, as discussed above. We note that out 
of 15 galaxies plotted in figure[8l only two (with Sersic indices 2.3 and 
2.4) would not have entered the "early-type sample" if using a Sersic 
index threshold of 2.5. 



termination is systematically offset (possibly due t o an offset 
with respect to the stellar mass estimates used in IShen et all 
(2003), rather than because of environ mental effects, see discus- 
sion in the original paper). Using the IValentinuzzi et al.l (1201 Ol) 
relation as the local reference thus reduces the difference be- 
tween our sample and local massive (cluster) early-types (me- 
dian r/r-=o=0.64±0.08). 

In addition, Marasto^ (l2005h has pointed out how a differ- 
ent treatment of the thermally pulsing asymptotic giant branch 
(TP-AGB) phase of stellar evolution may introduce additional 
(to those mentioned in section 14.2b systematics in the stellar 
masses, whose significance depends on the age of the stellar 
populations. This may be particularly relevant when compar- 
ing early-type galaxies at different redshifts, because while local 
early-types are too old for TP-AGB phase making a significant 
difference, at higher redshifts they inevitably become closer to 
hosting the intermediate age stellar populations for which a dif- 
ferent TP-AGB phase treatment may have an important effect. 
For the sample relev ant to Figure|8] stellar masses estimated with 
the MarastonI (12005 ) models are lower by about 0. Idex than the 
masses plotted, estimated with Bruzual & Chariot (2003). Using 
these lower masses, the size of galaxies in Figure |8] are ~ 50% 
and ~ 70% of those of e arly-types at z^O, based on t he lo- 
cal relations bv lShen etaTI {2003) and IValentinuzzi et al.l (1201 Oh 
(r/r.=o=0.5±0.1 and 0.72±0.1, respectively). We note that if we 
use the updated version of the Bruzual & Chariot ( 2003) models 
(often referred to as CB07), incorporating an improved prescrip- 
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Fig. 8. The stellar mass vs size relation for massive red- 
sequence early-type galaxies in the core of XMMU J2235. 
Filled symbols are spectroscopic members, empty symbols are 
sources with photo-z within 3 c r of the clus t er red shift. The black 
line and shaded area show the IShen et aT ] (120031) determination 
of the local stellar mass vs size relation and its scatter (Icr), 
bas ed on a sainple of SDSS early-types. The dotted line shows 
the Shen et aP (l2003h relation shifted by a factor 2 in size. The 
gray line with errorbars shows the loc al mass vs size relation 
determined bv IValentinuzzi et"ail (120101) for a sample of nearby 
cluster galaxies. 
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tion for the treatment of the TP-AGB phase, stellar masses for 
the Figure [8] sample are lower (in median) by 0.03dex; CB07 
mass estimates for this sa mple range from almost 3 0% lower to 
-10% higher compared to lBruzual & CharlotI (l2003i) . 

A combination of systematics in stellar masses and mismatch 
with the local reference relation might thus shift the (median) 
difference in size of XMMU J2235 and local early-types be- 
tween a factor slightly over 2 and a factor 1 .4. 

We note that other effects might also systematically bias 
the stellar masses and sizes used above. Again, we recall 
that our sizes were estimated in the restframe U band, while 
we are comparing with l ocal galaxy sizes m easured in the z 
dShen et al. 2003) and V (IValentinuzzi et al! 2010.) passbands. 



Based on previous work (jMclntos h et al. 2005 : Barden et al.l 



120051; iTruiillo et al]|2007l; ISzomoru et al.ll2010l) . correcting our 
sizes to the restframe z/V bands might reduce our sizes by 
a factor of up to 30%, and thus increase the median differ- 
enc e in size with local count erparts to up to a fact or ~ 3 (us- 
ing Bruzual & Charloj (l2003h stellar masses and the IShen et aTl 
local reference). 
We also remind the reader that the stellar masses used here 
were derived assuming a fixed, solar metallicity. Should the 
metallicity of the galaxies in Figure [8]be lower than solar, their 
stellar masses would be higher, by a factor about 15% and 2% 
for metallicites of 0.004 and 0.008, respectively. If their metallic- 
ity is supersolar, their stellar masses would be lower, by a factor 
-30% for a metallicity of 0.05. 

Finally, we note that the sizes we used were estimated based 
on a Sersic profile fit with varying nsersic- Therefore, for the sam- 
ple of early types plotted in Figure [8] nsei sic is not fixed to 4, but 
ranges between ~ 2.5 - 6 with an average nsersic ~ 3.5. Because 
of the correlation between the estimated ngersic and size, an ob- 
ject with a nsersic^'^ profile which is fitted with a lower Sersic in- 
dex (possibly due to the faintness of the source) will also have its 
size biased to a lower value. In order to estimate the relevance of 
such an effect, we fitted all sources in Figure|8]with nsersu- fixed 
to 4. This produces a very mild difference (on average <15%) in 
the estimated sizes, and overall negligible difference (<5%) in 
the evolution factors quoted above. 

Although the above results imply an evolution of the stellar 
mass vs size relation by a factor of about 1 .4 to more than 2, this 
evolution does not necessarily imply that the invidual galaxies 
in our early-type sample will evolve by this factor by redshift 
zero. Several studies showed how the population of early-type 
galaxies is increased, as a function of cosmic time, by the evolu- 
tion of late-type, star-forming galaxies into quiescent early-types 
(e.g., among many other s. Bell et al. 2004; Pannella et al. 2006; 



(e.g., among man y 
iFaber et al.ll2007H l 



iBrown et al.1120071; IFranx et al.ll2008L and ref- 



erences therein). It is thus very likely that a significant fraction 
of early-types at z=0 were still forming stars at z= 1 .4 (or did not 
yet have an early-type morphology), and were added to the lo- 
cal "reference" early-type sample at later cosmic times. Together 
with the present-day correlation (at fixed mass) between galaxy 
size and age of its host stellar populations, this may indeed 
suggest that the comparison of high-redshift early-type samples 
with local samples might artificially produce a "size evolution" 
signature, because the high-redshift sample does not account 
for the whole pr ogenitor populatio n of the local sample (se e 
discussi on in e.g. l Franx et al. 1 (|2008l); Ivan der Wei et all (l2009h: 



[Valentin uzzi et al.r(l2010h;lB ernardi et al. (2010h; IWimams etalj 
(12010)). According to van der Wei et al.1 (12009) (their Figure 5) 
and Valentinuzzi et al. (2010) (their Figure 4), the systematic 
exclusion from the high-redshift sample of the progenitors of 
younger local early-types, which tend to have larger sizes, might 



produce a difference in size of a factor as large as 1.3 at the red- 
shift of XMMU J2235. 

A definite conclusion on the actual relevance of size evolu- 
tion for our sample of massive cluster early-types is thus pre- 
cluded by the possible biases mentioned above. It remains clear 
that the significant range of galaxy sizes at all redshifts requires 
much larger samples in order to draw a conclusive picture of the 
mass and environmental dependence of galaxy size evolution. 

Finally, for the sake of completeness, we note that by com- 
paring the sizes of our spectroscopic members of late-type mor- 
phology w ith the typi c al siz e of late-type galaxies in the nearby 
Universe dShen et alJ ((2003), as well as iMaltbv et al] (120101) 
specifically in "cluster" and "cluster core" environments), we 
cannot find evidence of significant size evolution at fixed stel- 
lar mass for disk-dominated galaxies. The formal difference we 
find by comparing our late-type s izes at masses 10'° - IO^'Mq 
with cluster galaxies at z~0.2 (Maltbvetal. 2010) would be 
up to ~15% (for our standard Bruzual & Chariot (2003) stellar 
masses), but given the large uncertainties inherent in the com- 
parison (as discussed above), and the fact that our spectroscopic 
late-type sample suitable for this comparison is small (< 10 
galaxies) and incomplete, we cannot derive any definite conclu- 
sion. 



6. The stellar mass function 

Since the observed Kj-band light probes the rest-frame z band of 
galaxies at the cluster redshift, the K^-band LF derived in section 
[3]can be used to estimate the stellar mass function of galaxies in 
the central region (r < rsoo) of XMMU J2235. 

We note that, while a galaxy NIR luminosity has a milder 
sensitivity to recent star formation and dust attenuation as com- 
pared to luminosity at shorter wavelengths, stellar masses de- 
rived from the NIR luminosity alone may still suffer from sys- 
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Fig. 9. The galaxy stellar mass function in the central region 
of XMMU J2235, as derived from the Kj-band LF (soUd sym- 
bols). Empty black circles show the stellar mass function de- 
rived from individual SED-estimated stellar masses of spec- 
troscopic and photo-z retained members (se e text for details). 
Gray symbols show previou s determinations dKodama & BoweJ 
l2003blStrazzullo et al.ll2006l) at z ~ 1 (squares) and z ~ 1.2 (cir- 
cles). All mass functions are arbitrarily rescaled. The dotted line 
shows the mass completeness limit as estimated from the K, im- 
age lOcr completeness. 
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tematics of up to a factor ~2 depending on different mass-to- 
light (M/L) ratios of different galaxies. Based on stellar masses 
of sources with photo-zs within 2cr of the cluster redshift, we 
calibrated the stellar mass vs observed Kj-band luminosity re- 
lation finding a scatter of ~ 0.3dex for the whole sample and 
of ~ O.lSdex when considering blue and red galaxies separately 
(using a color threshold of z-H^l.5). 

We therefore applied a statistical correction for different M/L 
ratios of cluster galaxies in the Kj-band selected sample, and de- 
rived the mass function and its uncertainties as follows. Firstly, 
we estimated the fraction of red and blue cluster galaxies as a 
function of K^-band luminosity based on the 2cr photo-z retained 
sample, based on both the z-H and z-J colors, yielding consistent 
results. Because of the ~100% completeness and ~5Q% contam- 
ination of this sample, we may expect that the estimated rela- 
tive contributions of red and blue cluster galaxies are possibly 
biased toward a higher blue fraction, since the interlopers con- 
taminating our sample will be, on average, preferentially bluer 
than cluster galaxies. This would result in a possibly lower M/L 
ratio used to statistically convert the K^-band light to stellar mass 
(and thus lower masses). 

Once the relative contributions of red and blue galaxies, and 
the typical M/L ratios of these two populations were determined, 
the K-band LF was translated into a stellar mass function by 
means of 1000 realizations taking into account the error on the 
LF (as plotted in Figure[3]l, the estimated scatter in the M/L ratios 
of red and blue galaxies, and the error in the red/blue galaxy 
fractions. 

The derived stellar mass function is shown in Figure|9] Solid 
black symbols and errorbars show the median and 16-84f/! per- 
centiles of the distribution obtained through the 1000 realiza- 
tions. For comparison, we also show (empty black symbols) the 
stellar mass function that is directly obtained from the individual 
SED-estimated stellar masses of the spectroscopic members plus 
photo-z retained sources randomly sampled assuming a ~50% 
contamination by interlopers. These two largely independent de- 
terminations of the stellar mass function in the central region of 
XMMU J2235 are perfectly consistent. The shape of the mass 
function determined here is in very good agreement with (ar- 
bitrarily res caled) previous measurements in clusters at z ~ 1 
and z ~ 1 .2 ( Kodama & Bower 2003:'Strazzull o et al.l200 ^, both 
scaled to a Kroupa IMF), once again suggesting an early assem- 
bly of massive galaxies at least in overdense environments, with 
the high-mass end population essentially in place at one third of 
the Hubble time. 

The area considered here lies within r <90"=765kpc, 
which is very close t o the estimated r ^inn for XMMU 
J2235 (rsoo ~0.75Mpc, iRosati et all l2009h . The total pro- 
jected stellar mass within this area, inferred from integrat- 
ing the Kj band LF and converting the total K-band luminos- 
ity (LK,5oo,pmj-^''L*T(2+a)-ll+2 IO'^Lq) to stellar mass, is 
M .5fflr.5,500,praj— 6^^ ■ lO'^M©, with errors including the uncertain- 
ties on the LF parameters and the scatter in the galaxy M/L ra- 
tios. For comparison, the stellar mass in galaxies more massive 
than ~ 1O'"M0, statistically determined from the photo-z re- 
tained sample within the same area (see mass function in Figure 
HI), is ~ 3 ■ lO^^Mo^ 

Following iRosati et"al] ( l2009l) the total projected mass 
within this radius is M,ot^5oo,proj = 7 + 2 ■ lO'^^M©. The ratio of 
total mass to K-band luminosity in the projected area within rsoo 
is thus M,o,,5oo,pro;/L;r,500,proj=60±20 Mq/Lq, while the stellar 
mass fraction is Ms„„,50(),/,rayM,o,_5oo,proj=0.009+[J{j|^7, in close 
agreement with determinations for massive clusters at lower 



redsh ifts (e.g.. i Lin et al.ll200l lEttori et alJl2009t iGiodini et al.1 
120091: lAndreonl2010h . 

7. Summary 

We presented a study of galaxy populations in the core of the 
massive cluster XMMU J2235 at z=1.39, making use of high 
quality VLT and HST multi-wavelength photometry in seven 
passbands from U to Kj, sampling the rest-frame SED of cluster 
members from the NUV to the NIR. 

We derived luminosity functions in the z, H, and K^ bands, 
approximately corresponding to the rest-frame U, R, and z band. 
These extend the study of the galaxy luminosity function in mas- 
sive X-ray luminous clusters beyond z ~ 1.3, and are among the 
deepest determinations at z > 1 . All three determinations of the 
luminosity function, probing down to ~M*-i-2 or M*+4 depend- 
ing on the passband, are consistent with having a flat faint-end 
slope and a characteristic magnitude M* close to passive evolu- 
tion predictions of the M* of local massive clusters with a for- 
mation redshift z > 2. 

The color-magnitude and color-mass diagrams in the core 
region of XMMU J2235 show evidence of a tight (intrinsic scat- 
ter <0.08) red sequence of massive galaxies, with overall old 
stellar populations (z/ > 2), no evidence of significant on-going 
star formation (SFR< 0.2MQ/yr), and generally early-type mor- 
phology. The spectroscopic sample is essentially complete at 
the high-mass end, and strongly suggests that the most massive 
cluster galaxy populations in the core of this cluster are already 
dominated by early-type galaxies, both in terms of galaxy struc- 
ture and of overall spectro-photometric properties. 

Star formation appears effectively quenched at masses higher 
than ~ 6- 1O'"M0. Also, active star formation is suppressed in the 
very central regions, with all spectroscopically confirmed star 
forming cluster members located at r > 250kpc from the BCG. 

In agreement with previous work at lower redshifts, these 
data point toward an early assembly of massive cluster galaxies, 
not only in terms of the formation of their stars, as suggested by 
their broad-band SEDs, as well as by their spectra (Rosa ti et al.l 
2009), but also of the assembly of their stellar mass, as suggested 
by the SED-derived stellar masses of individual sources, by the 
close-to-passive evolution of the luminosity function bright- 
end, and by the galaxy stellar mass function. As an additional 
evidence that XMMU J2235 is in a very advanced evolutionary 
stage already at this redshift, we estimate a stellar mass fraction 
within rsoo of about 1%, similar to local massive clusters. 

On the other hand, comparing the size of these massive red- 
sequence galaxies to the size of similarly massive local early- 
types, suggests a possible size evolution of up to a factor ~ 2. 
This implies that, in spite of the overall early assembly of these 
sources, room is left for processes like minor (and likely dry) 
merging to shape the structural properties of these objects to re- 
semble those of their local counterparts, without substantially 
affecting their stellar mass or galaxy populations. We discussed 
how the role of possible systematics makes it still difficult to 
draw firm conclusions on the magnitude of such a size evolution, 
and on whether this evolution might be different from the one 
estimated in field galaxy samples over the same redshift range. 
Much larger samples of high-redshift (cluster) early-type galax- 
ies are needed to derive a conclusive picture of the mass and 
environmental dependence of their size evolution. 

Future observations at IR, mm and radio wavelengths, prob- 
ing dust-unbiased star formation and the cold gas component, 
will enhance our understanding of the mass growth, star forma- 
tion, and structural evolution of high-redshift cluster galaxies. 
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